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UPDATE  TO  THE  FREE  FLOATING  TWO-DIMENSIONAL  EXTENSIBLE  CABLE  SYSTEM 

MODEL  (FF2E) 


PURPOSE 


The  purpose  of  this  report  is  to  explain  the  modifications  made  to  the  FF2E 
computer  program  to  more  accurately  predict  operational  behavior  of  a  drifting 
sonobuoy.  Each  modification  will  be  discussed  in  detail  including  the  new 
data  inputs  required  to  execute  the  program. 

BACKGROUND 


The  FF2E  computer  program,  references  (1)  and  (2),  was  developed  to 
predict  the  steady- state  behavior  of  a  cable  system  in  an  ocean  environment. 
The  program,  as  referenced  in  sonobuoy  specifications,  is  used  to  predict  the 
operating  conditions  of  a  sonobuoy.  The  two  conditions  generally  emphasized 
are  the  operating  depth  of  the  hydrophone  and  the  velocity  of  the  hydrophone, 
relative  to  the  current  at  the  hydrophone's  operating  depth.  For  the  purpose 
of  this  report,  the  relative  hydrophone  velocity  will  be  referred  to  as  the 
"flow  velocity".  The  flow  velocity  for  a  sonobuoy  design,  in  a  given 
two-dimensional  current  profile,  can  be  predicted  as  a  function  of  the  drag 
characteristics  of  the  sonobuoy  system.  The  predicted  flow  velocity  is  used 
to  establish  the  test  scenario  for  evaluating  mechanically-induced  noise.  The 
level  of  mechanically-induced  noise  is  a  function  of  the  flow  velocity; 
therefore,  an  accurate  flow  prediction  is  required.  With  increasing  emphasis 
being  placed  on  reducing  mechanically-induced  noise,  a  more  accurate  model  of 
a  drifting  sonobuoy  is  required. 

OBJECTIVE 


The  objective  of  the  task  was  to  develop  a  more  accurate  means  of  modeling 
the  behavior  of  a  drifting  sonobuoy.  The  task  was  divided  into  four  major 
areas: 


a.  Develop  a  subroutine  which  would  find  the  true  drag  force  on  a  surface 
unit  given  the  flow  velocity  past  the  surface  unit. 

b.  Incorporate  a  change  to  the  program  which  will  allow  the  use  of  the 
analytical  formula  for  the  90%  ocean  current  profile  of  reference  (3)  to 
calculate  flow  velocity  at  any  given  depth. 

c.  Add  a  subroutine  which  would  compensate  for  the  non-linear  stretch 
characteristics  of  a  cable  or  compliance  in  the  sonobuoy  suspension  system. 

d.  Make  a  modification  to  the  subroutine  "BODY"  which  corrects  an 
inherent  interpolation  error. 
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SUMMARY  OF  CHANGES 


a.  An  Interpolation  scheme  has  been  Included  in  the  subroutine  "VAR" 
which  allows  the  true  drag  force  on  a  surface  unit  to  be  found.  This  is 
accomplished  by  interpolation  of  empirically  measured  values  relating  drag 
forces  relative  to  flow  velocities  on  the  surface  unit.  If  no  table  of  values 
Is  to  be  input,  the  program  will  operate  as  previously  written. 

b.  A  modification  to  the  subroutine  “CUR"  has  been  incorporated  which 
will  allow  direct  calculation  of  flow  velocity  given  any  depth  in  the  90% 
ocean  current  profile.  When  this  profile  is  used,  input  of  current  velocity 
data  points  is  not  required.  If  the  90%  profile  is  not  used,  the  program  will 
operate  as  previously  written,  by  Interpolation  between  current  velocity  data 
points. 


c.  The  subroutine  "STRETH",  identified  by  reference  (4),  has  been  added  to 
the  program  to  account  for  the  non-linear  stretching  characteristics  of  cables 
and  compliances.  This  subroutine  utilizes  an  input  table  of  tension  values 
versus  AE  values,  where  "A"  is  equal  to  the  cross-sectional  area  and  "E"  is 
equal  to  Young's  Modulus.  If  a  cable  has  a  constant  value  of  AE  the  program 
will  operate  as  previously  written.  The  input  for  AE  has  been  altered  to 
allow  for  the  input  of  tables. 

d.  A  modification  has  been  made  to  the  subroutine  "BODY"  which  corrects 
an  error  in  output  when  a  velocity  or  angle  value  falls  within  the  first 
segment  of  the  lift  drag  table. 

DISCUSSION 


Because  of  Fleet  operational  requirements,  increased  emphasis  has  been 
placed  on  reducing  mechanically  induced  noise  in  sonobuoys.  Therefore, 
emphasis  is  being  placed  on  upgrading  and  refining  the  analytical  tools  used 
to  evaluate  the  suspension  systems. 

Due  to  this  fact,  modifications  to  the  FF2E  program  have  been  incorporated 
to  create  a  more  accurate  prediction  of  a  drifting  sonobuoy.  In  the  following 
paragraphs  each  modification  will  be  discussed.  Included  in  this  report  is  an 
improved  data  input  scenario  which  contains  the  added  input  parameters  for  the 
modifications.  In  Appendix  I  there  is  a  complete  listing  of  the  improved  FF2E 
program. 

SUBROUTINE  "VAR" 

In  the  present  version  of  the  FF2E  program,  drag  of  the  surface  unit  is 
determined  using  the  flow  velocity  past  the  unit  at  a  depth  of  one-half  the 
draft,  where  the  draft  is  equal  to  the  distance  from  sea  surface  to  the  bottom 
of  the  surface  unit.  As  the  program  iterates,  it  searches  for  an  equilibrium 
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of  the  buoyancy  of  the  upper  unit  and  the  effective  weight  of  the  hydrophone 
and  suspension  system  by  varying  the  drift  velocity  and  sonobuoy  draft. 

Because  the  FF2E  program  uses  a  mathematical  equation  to  model  the  drag  force 
of  the  upper  unit,  the  calculated  drag  force  may  not  exactly  model  the 
empirically  derived  drag  force  of  the  actual  hardware.  This  inherent  error 
creates  a  false  equilibrium  of  forces  and  ultimately  could  effect  the  flow 
velocity  of  the  hydrophone. 

The  subroutine  "VAR"  was  added  to  allow  the  true  drag  of  the  surface  unit 
to  be  found  given  the  relative  flow  of  the  upper  unit.  This  subroutine 
utilizes  a  reference  table  of  values  which  contain  measured  drag  forces  versus 
flow  velocities.  In  each  iteration  of  the  subroutine  “STEADY"  the  subroutine 
"VAR"  is  called  with  values  of  drag,  flow,  draft,  diameter,  and  drag 
coefficient.  Using  the  value  of  flow  equal  to  the  flow  velocity  at  one-half 
the  draft,  the  subroutine  finds  the  drag  force  by  interpolating  between 
empirically  derived  data  points.  This  value  of  drag  force  is  then  compared 
with  the  calculated  drag  force.  If  the  difference  between  these  two  values  is 
less  than  .001  lbs.,  the  program  will  return  and  continue  with  the  iterations 
of  the  subroutine  "STEADY".  For  differences  greater  than  .001  lbs.,  the 
subroutine  calculates  the  surface  unit  drag  coefficient  based  on  the 
empirically  measured  drag  force  value.  This  coefficient  is  then  used  to  vary 
the  original  surface  unit  drag  coefficient  until  correlation  between  the 
empirically  measured  drag  force  and  the  calculated  drag  force  is  achieved. 

DATA  INPUT  CHANGES  FOR  SUBROUTINE  "VAR" 

Three  variables  have  been  added  to  the  data  input  scheme.  The  first  is 
NFOSB.  NFOSB  is  the  number  of  drag  forces  to  be  entered  in  the  drag  versus 
flow  velocity  table.  A  maximum  of  15  values  may  be  entered  in  the  table.  The 
second  variable  is  FOSB(k),  where  k=l  to  NFOSB.  FOSB(k)  corresponds  to  the 
drag  force  values  to  be  input  in  the  table.  FOSB(k)  should  range  from  0  lbs. 
to  a  magnitude  which  is  greater  than  the  expected  drag  on  the  surface  unit. 

The  final  variable  is  VOSB(k),  which  is  the  flow  velocity  in  knots  that 
corresponds  to  the  force  values  entered  in  FOSB(k).  If  no  table  is  to  be 
entered  for  the  surface  unit  drag,  then  NFOSB  should  be  set  equal  to  zero,  and 
no  values  for  FOSB(k)  and  VOSB(k)  are  entered. 

DATA  OUTPUT  FOR  SUBROUTINE  "VAR" 

The  program  will  print  out  the  corresponding  value  of  FOSB  and  VOSB  in  a 
table  format.  After  the  program  has  reached  the  steady  state  solution,  it 
will  print  out  the  drag  force  on  the  surface  unit  as,  “the  computed  horizontal 
component  of  tension".  Also,  the  program  prints  out,  "the  CD  of  the  surface 
unit",  and  "the  flow  past  the  surface  unit". 
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SUBROUTINE  CUR 


In  the  present  FF2E  program,  subroutine  “CUR"  is  called  to  perform  a 
linear  interpolation  of  the  input  current  profile  points.  The  purpose  of  this 
subroutine  is  to  find  a  flow  velocity  for  a  given  depth  in  the  ocean  profile. 
Oepending  on  the  step-size  and  number  of  current  profile  points  entered,  the 
amount  of  error  in  the  interpolation  of  depths  close  to  the  surface  will  vary, 
since  the  current  profile  is  not  a  straight  line  at  shallow  depths. 

In  the  modified  program,  the  analytical  formula  which  describes  the  90% 
profile,  has  been  added  to  subroutine  "CUR“.  If  the  user  wishes  to  use  the 
90%  profile,  then  the  number  of  current  profile  points,  NCUR,  must  be  set 
equal  to  zero  and  XX ( I )  and  Y Y { I )  values  are  not  entered.  If  NCUR  is  not 
zero,  then  the  values  for  the  current  profile  must  be  entered. 

SUBROUTINE  "STREW 

The  present  FF2E  computer  program  assumes  that  the  cable  and  compliance  of 
a  sonobuoy  suspension  system  will  have  linearly  stretching  characteristics.  In 
the  past,  the  problem  of  a  non-linear  stretching  cable  has  been  over-ridden  by 
inputting  the  stretched  length  of  the  cable  into  FLC(k)  and  choosing  an  AE 
value  of  high  magnitude. 

In  the  modified  program,  subroutine  "STRETH",  identified  by  refrence  (4), 
has  been  added  to  account  for  cables  with  non-linear  stretch  characteri sties. 
The  subroutine  utilizes  a  reference  table  of  values  which  relates  cable 
tension  to  a  value  of  AE.  The  AE  value  is  calculated  from  the  following 
equation: 


AE  =  Lo  (T-Tref) 


L-Lo 


Where,  2 

A  =  cross  sectional  area  (ft*) 

E  =  Young's  Modulus  (lb/ft  ) 

L  =  stretched  length  of  cable  (ft) 

Lo  =  unstretched  length  of  cable  (ft) 

T  =  tension  applied  to  cable  (lb) 

Tref  =  reference  tension  of  unstretched  cable  (lb) 


DATA  INPUT  CHANGES  FOR  "STRETH" 

Four  variables  have  been  added  to  the  input  data  parameters.  The  first  is 
NAE  which  corresponds  to  the  number  of  points  to  be  entered  for  the  relation 
of  tension  to  AE  for  each  cable  and  compliance.  NAE  has  a  range  from  1  to  15. 
FAE(I,K)  and  AE( I ,K )  are  the  corresponding  tension  and  AE  values  for  a  cable 
over  which  the  tension-strain  curve  is  non-linear,  where  1=1  to  NCAB  and  K=1 
to  NAE.  For  cables  with  a  linear  tension-strain  curve,  NAE  will  be  set  equal 

to  1,  and  the  fourth  variable,  AE ( I , 1 )  would  be  equal  to  a  constant  AE  value. 
Values  for  AE  should  be  entered  for  the  entire  range  of  tensions  over  which 
the  tension-strain  curve  is  non-linear.  For  tension  values  outside  the  range 
of  the  table,  AE  will  be  set  equal  to  the  endpoint  value. 
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DATA  OUTPUT  CHANGES  FOR  "STREW 

The  program  will  print  out  the  corresponding  Tension-AE  values  in  tabular 
format  as  they  are  entered  in  the  input.  For  cables  with  linear  stretch 
characteristics  the  constant  value  of  AE  will  be  printed  out. 

SUBROUTINE  "BODY" 

A  modification  to  the  subroutine  “BODY",  reference  (4),  has  been  provided 
which  corrects  an  error  in  this  subroutine.  Previously,  when  a  velocity  or 
angle  value  fell  within  the  first  segment  of  the  lift  drag  table,  the  first 
value  of  the  table  was  chosen  rather  than  interpolating  between  the  first  and 
second  values  of  the  table.  As  an  example,  assume  the  following  table  is  in 
the  program: 


ANGLE 


0 

10 

20 
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L 

D 

L 

D 

L 

D 

E 

1 

- 

4 

1 

7 

4 

10 

7 
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0 

2 

- 

5 

2 
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5 

11 

8 
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- 

I 

3 

- 

6 

3 

9 

6 

12 

9 

T 

Y 

The  new  data  points  may  be  computed  as  follows: 


STEP  1:  Determine  lift  at  the  desired  velocity  and  at  a  given  angle  by 
Interpolation 

(V3  '  V  (V2  '  V  (V2  "  V  +  Li  =  12 


where  -  desired  velocity 

=  next  lowest  velocity  data  point 
=  next  highest  velocity  data  point 
L?  =  calculated  lift  for  desired  velocity 
if  =  lift  value  for  next  lowest  velocity 
Lj  =  lift  value  for  next  highest  velocity 


Example: 


(2-1)  (1.5  -  1)  +  4  =  L0  =  4.5  for  angle  =  0°,  and  velocity  =  1.5 
- -  2 

(2  -  1 )  (1.5  -  1 )  +  7  =  l "2  *  7.5  for  angle  =  10°,  and  velocity  =  1. 


(8  -  7) 
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STEP  2:  Determine  lift  for  the  desired  angle  at  the  desired  velocity  by 
Interpolation: 

(L'2  -  L2)  (A2  -  Aj)  +  L2  =  L 


where  A.  =  next  lowest  angle  in  table 
Ai  =  next  highest  angle  in  table 

=  value  for  lift  at  the  desired  velocity  and  angle 

Example:  (7.5  -  4.5)  (10  -  0)  +  4.5  =  6  for  angle  =  5°  and  velocity  =  1.5 

5-0 

STEP  3:  Compute  drag  by  substituting  "D"  for  "L"  and  repeating  steps  1  and  2. 

The  new  data  points  for  the  example  are  found  to  be  lift  =  6  and  drag  =  3 
for  velocity  =  1.5  and  angle  =  5°. 


6 


NAC-TR-2359 


INPUT  DATA  FORMAT  AND  DESCRIPTION 


A  description  of  the  input  parameters  of  the  modified  FF2E  program  is 
shown  below: 

CARD  1 

COLUMN  1-3 

1.  NCASES  -  number  of  sets  of  data  to  be  ran.  NCASES  is  set  equal  to  1  for 
one  set  of  data  to  be  ran.  NCASES  is  an  integer. 

CARD  2 

COLUMN  1-3 

2.  NCUR  -  number  of  current  profile  points  to  be  entered.  The  maximum  number 
of  points  is  30.  If  NCUR=0,  the  90%  profile  is  assumed.  NCUR  is  an  integer. 

COLUMN  4-6 

3.  NCAB  -  number  of  cable  segments  bounded  by  endpoints  or  system  components. 
NCAB  is  an  integer. 

COLUMN  7-9 

4.  NHPHS  -  number  of  hydrophones.  NHPHS  is  an  integer. 

COLUMN  10-12 

5.  NTAB  -  number  of  lift/drag  tables.  If  no  tables  are  used  then  NTAB=0. 
NTAB  is  an  integer. 

CARD  3 

COLUMN  1-12 

6.  DAI  -  diameter  of  buoy  antenna  in  inches.  DAI  is  a  floating  point  number. 


COLUMN  13-24 


7.  LA  -  length  of  buoy  antenna  in  feet.  LA  is  a  floating  point  number. 
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COLUMN  25-36 

8.  CDA  -  drag  coefficient  in  air  of  buoy  antenna.  CDA  is  a  floating  point 
number. 

COLUMN  37-48 

9.  T8H  -  horizontal  component  of  tension  in  pounds  applied  to  lower  unit. 

TBH  is  a  floating  point  number. 

10.  TBV  -  vertical  component  of  tension  is  pounds  applied  to  lower  unit.  TBV 
is  a  floating  point  number.  TBH,  TBV  should  equal  0.  for  free  floating  system. 

CARD  4 

COLUMN  1-12 

11.  DB  -  diameter  of  cylindrical  buoy  in  feet.  DB  is  a  floating  point  number. 
COLUMN  13-24 

12.  LB  -  length  of  cylindrical  buoy  in  feet.  LB  is  a  floating  point  number. 
COLUMN  25-36 

13.  CDB1  -  drag  coefficient  of  submerged  portion  of  surface  buoy.  CDB1  is  a 
floating  point  number. 

COLUMN  37-48 

14.  CDB2  -  drag  coefficient  of  surface  buoy  above  sea  surface.  CDB2  is  a 
floating  point  number. 

COLUMN  49-60 

15.  WB  -  weight  of  surface  buoy  in  air  in  pounds.  WB  is  a  floating  point 
number. 

COLUMN  61-72 

16.  UWINDK  -  velocity  of  wind  in  knots.  UWINDK  is  a  floating  point  number. 

CARO  5 

COLUMN  1-3 

17.  NFOSB  -  number  of  values  on  the  force  vs.  flow  curve.  NFOSB  =  0  to  15 
and  is  an  integer. 
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CARD  6 

COLUMN  1-12,  13-24,  etc. 

18.  FOSB(K)  -  force  values  on  the  force/flow  curve  for  the  surface  unit  in 
pounds.  FOSB(k)  are  floating  point  numbers,  where  K=l,  NFOjB. 

CARD  7 

COLUMN  1-12,  13-24,  etc. 

19.  VOSB(K)  -  flow  values  on  the  force/flow  curve  for  the  surface  unit  in 
knots.  VOSB(K)  are  floating  point  numbers,  where  K=l,  NFOSB. 

Note:  For  variables  18  and  19,  a  maximum  of  6  values  may  be  read  in  on 
one  card.  As  many  additional  cards  as  required  may  be  used. 

CARD  8 

COLUMN  1-12 

20.  CDAPK  -  drag  coefficient  times  the  normal  drag  area  of  package  hung  just 
below  the  surface  buoy  in  feet  squared.  CDAPK  is  a  floating  point  number. 

COLUMN  13-24 


21.  WPAK  -  wet  weight  of  package  hung  just  below  the  surface  buoy  in  pounds. 
WPAK  is  a  floating  point  number. 

CARD  9 

COLUMN  1-12,  13-24,  etc. 

22.  FLC(K)  -  length  of  each  cable  in  feet.  FLC(K)  is  a  floating  point 
number,  where  K=l,  NCAB. 

CARD  10 

COLUMN  1-12,  13-24,  etc. 

23.  NPR(K)  -  number  of  segments  each  cable  is  divided.  NPR(K)  are  integer 
values,  where  K=l,  NCAB. 
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CARO  11 

COLUMN  1-12,  13-24,  etc. 

24.  DCI(K)  -  diameter  of  each  cable  In  inches.  DCI(K)  are  floating  point 
numbers,  where  K=l,  NCAB. 

CARD  12 

COLUMN  1-12,  13-24,  etc. 


25.  CVFAC(K)  -  conversion  factor  for  non-circular  cross  section  of  cable. 
CVFAC  is  a  floating  point  number,  where  K=l,  NCAB. 

CARD  13 

COLUMN  1-12,  13-24,  etc. 

26.  WC(K)  -  weight  per  unit  length  of  each  cable  in  water  in  pounds  per  foot. 
WC(K)  are  floating  point  numbers  where  K=l,  NCAB. 

CARD  14 

COLUMN  1-12,  13-24,  etc. 

27.  CDC(K)  -  coefficient  of  drag  on  each  cable  normal  to  the  cable.  CDC(K) 
are  floating  point  numbers,  where  K=l,  NCAB. 

CARD  15 

COLUMN  1-12,  13-24,  etc. 

28.  TREF(K)  -  reference  tension  in  each  cable  in  pounds.  If  FLC(K)  is  an 
unstretched  length,  set  TREF(K)  =  0.  TREF(K)  are  floating  point  numbers. 

CARD  16 

COLUMN  1-12,  13-24,  etc. 

29.  (K)  -  Poisson's  ratio  for  each  cable.  P(K)  are  floating  point  numbers. 
CARD  17 

COLUMN  1-3 

30.  NAE(K)  -  number  of  points  on  the  tension/AE  curve  for  each  cable.  NAE(K) 
are  integers.  For  a  cable  without  a  table  NAE=1. 


V. 

i 
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CARD  18 

COLUMN  1-12,  13-24,  etc. 

31.  AE(K,1)  -  AE  value  for  cable  with  linear  stretch  characteristics  with 
NAE=1.  AE(K,1)  is  a  floating  point  number. 

CARD  18a 

COLUMN  1-12,  13-24,  etc. 

32.  FAE(K,I)  -  tension  values  over  which  AE  is  changing  in  pounds.  FAE(K,I) 
are  floating  point  numbers,  where  K=l,  NCAB  and  1=1,  NAE. 

CARD  18b 

COLUMN  1-12,  13-24,  etc. 

33.  AE(K,I)  -  corresponding  values  of  AE  relative  to  the  values  of  card  18a. 
AE(K,I)  are  floating  point  numbers. 

CARD  19 

COLUMN  1-12,  13-24,  etc. 

34.  CDABD(K)  -  Drag  coefficient  times  normal  drag  area  of  body  at  the  bottom 
of  each  cable  in  feet  squared.  CDABD(K)  are  floating  point  numbers,  where 
K=1 ,  NCAB. 

CARD  20 

COLUMN  1-12,  13-24,  etc. 

35.  WBD(K)  -  wet  weight  of  body  at  the  bottom  of  each  cable  in  pounds. 

WBD(K)  are  floating  point  numbers,  where  K=l,  NCAB. 

CARD  21 

COLUMN  1-3,  4-6,  etc. 

36.  NBOD(K)  -  lift/drag  table  number  which  applies  to  the  body  at  the  bottom 
of  each  cable.  If  no  body  table,  then  NB0D=0.  NBOD(K)  are  integer  numbers, 
where  K=l,  NCAB. 

CARD  22 

COLUMN  1-2,  13-24,  etc. 

37.  XX(I)  -  depth  of  each  current  point  in  feet.  XX ( I )  are  floating  point 
numbers,  where  1=1,  NCUR. 
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CARO  23 

COLUMN  1-12,  13-24,  etc. 

38.  YY(I)  -  current  velocity  at  each  point  corresponding  to  depths  of  XX ( I ) . 
Velocities  are  in  knots  and  are  floating  point  numbers,  where  1=1,  NCUR. 

Note:  If  NCUR=0,  then  values  for  XX ( I ) ,  Y Y ( I )  are  not  input. 

CARD  24 

COLUMN  1-3 

39.  NPHI(N)  -  number  of  angle  values  for  which  lift/drag  tables  are 
generated.  NPHI(N)  are  integer  numbers. 

COLUMN  4-6 

40.  NU(N)  -  number  of  velocity  values  for  which  lift/drag  tables  are 
generated.  NU(N)  are  integer  numbers. 

CARD  25 

41.  PHIM(N, I )  -  angle  in  degrees  of  suspension  cable  from  vertical  for  which 
lift/drag  forces  are  generated.  PHIM{N,I}  are  floating  point  numbers. 

CARD  26 

COLUMN  1-12,  13-24,  etc. 

42.  U(N,J)  -  Velocity  in  knots  for  which  lift/drag  forces  are  generated. 
U(N,J)  are  floating  point  numbers. 

CARD  27 

COLUMN  1-12,  13-24,  etc. 

43.  DDRAG(M)  -  drag  values  in  pounds  for  a  velocity  and  angle.  DDRAG(M)  are 
floating  point  numbers. 

CARD  28 

COLUMN  1-12,  13-24,  etc. 

44.  DLIFT(M)  -  lift  values  in  pounds  for  a  velocity  and  angle.  DLIFT(M)  are 
floating  point  numbers. 

Note:  For  variables  39  to  42,  N=1  to  NTAB,  1=1  to  NPHI(N),  and  J=1  to 
NU(N).  For  variables  43  and  44,  M=1  to  NP,  where  NP=NU(N)  times  NPHI(N). 
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APPENDIX  I 
FF2E  LISTING 
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